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ABSTRACT 

The HST ACS Coma Cluster Treasury Survey is a deep two passband imaging survey 
of the nearest very rich cluster of galaxies, covering a range of galaxy density envi- 
ronments. The imaging is complemented by a recent wide field redshift survey of the 
cluster conducted with Hectospec on the 6.5m MMT. Among the many scientific ap- 
plications for this data are the search for compact galaxies. In this paper, we present 
the discovery of seven compact (but quite luminous) stellar systems, ranging from 
■ M32-like galaxies down to ultra-compact dwarfs (UCDs) / dwarf to globular transition 

objects (DGTOs). 

We find that all seven compact galaxies require a two-component fit to their 
light profile and have measured velocity dispersions that exceed those expected for 
typical early-type galaxies at their luminosity. From our structural parameter analysis 
we conclude that three of the sample should be classified as compact ellipticals or 
M32-like galaxies, the remaining four being less extreme systems. The three compact 
ellipticals are all found to have old luminosity weighted ages (^12 Gyr), intermediate 
metallicities (-0.6 < [Fe/H] < -0.1) and high [Mg/Fe] (> 0.25). 

Our findings support a tidal stripping scenario as the formation mode of compact 
galaxies covering the luminosity range studied here. We speculate that at least two 
early- type morphologies may serve as the progenitor of compact galaxies in clusters. 

Key words: surveys — galaxies: clusters: individual: Coma — galaxies: dwarf — 
galaxies: structure 



1 INTRODUCTION 
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The term compact stellar system covers a range of possi- 
ble object types, characterised by small physical dimensions 
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and high surface brightness. The largest are the compact 
elhptical (cE) galaxies exemplified by M32. As a close com- 
panion of M31, its truncated radial profile has long been 
suspected to be due to tidal or other inte ractions with its 
giant neighbour ([Kind Il962l : lFabei]ll973l l. M32 itsel f has 
My — —16.5 mag and a half-light radius Re ^ 120 pc (|Kentl 
1 19871 ). Despite a luminosity in the dwarf galaxy regime, it 
has other properties more akin to those of a giant galaxy or 
a classical bulg:e, such as the presenc e of a massive central 
bla ck hole (|van der Marel et al.l [l997l ) and a high metallic- 
ity teo se et al.ll2005l ). reinforcing the suggestion that it has 
suffered major loss of mass. 

Few other galaxies with similar properties are known. 
The other usually quoted examples - somewhat brighter at 
- -18 mag, and with Re - 200-300 pc - are NGC 
4486B, a companion of M87 near the centre of the Virgo 
Cluster (Rood 1965), NGC 5846A (in th e NGC 5846 group ; 
lFaberlll973l ). two obie cts in Abell 1689 (iMieske et al.il2005l l 
and one in AbeU 496 (jChilingarian et al.ll2007l i^ 

Globular clusters (GCs) obviously also satisfy the re- 
quirements for compact stellar systems and are, of course, 
extremely numerous. Their typical luminosities (the turn 
over point in their luminosity functi on) are around My — 
—7.5 mag (|Di Criscienzo et"aDl2006l ) and the brightest ex- 
amples around our Galaxy [uj Cen) and M31 (Gl) reach 
My — —10 to —11 mag; the very populous GC systems 
around central cluster gala xies may exte nd slightly brighter, 
to -11.5 mag (Ha rris et al.ll2009l ). 

No compact object^] at luminosities between these two 
extreme s were known un t il the discovery in the Fornax 
Cluster (Hilker et al.' ll999l : IPrinkwater e t al. 2000) of what 
are generally referre d to as ultra-compact dwarfs (UCDs; 
IPhillipDS et alll200lh . These have My between -11.5 and 
— 13.5 mag and Re ~ 20-30 pc. Furt her similar exam- 
ples w ere found in the Virgo Cluster ([Prinkwater et all 



I2OO4I : IJones et al 



2006) and in a couple of large 



groups iFirth et al.1 12OO6I : lEvstigneeva et al.1 l2007al ) while 
iHasegan et al.1 (|2005 l added the slightly fainter, intermedi- 



ate objects which they call DGTOs, for d warf /globula r tran- 
sition objects (see also Drinkwater et al. 2005; Gre gg et al. 
I2OO9I) . iMarzke et all (j2006l ) and IChilingarian Mamon 
(|2008l ) have also suggested the existence, in Sloan Digital 
Sky Surve}Q(SDSS) data, of compact systems with lumi- 
nosities between UCDs and compact ellipticals. Between 
them, these new types of system now span essentially all of 
the luminosity space from normal globular clusters upwards. 
The majority appear to reside in reasonably dense environ- 
ments where interactions can be expected. Indeed even those 
found in relatively low density environments, such as the re- 
cently discovered UCD associated with the Sombrero galaxy 
fM104: lHau et al.ll2009l ) and M32 itself, have more massive 
companions. 

In the present paper we present the discovery of com- 
pact stellar systems in the very rich and dense Coma Clus- 
ter. The sample was observed as part of the Hubble Space 



^ We exclude blue compact dwarfs (BCDs) from this list, as they 
appear to be relatively diffuse underlyin g galaxies with very cora - 
pact central star forming regions, see e.g. lPapaderos et al.l (|2006l) . 
W e note, however, t hat so me BCDs, such as the object reported 
bv lTrentham et"al are analogous to blue cEs. 

^ http://www.sdss.org 



Telescope (HST) Treasury Survey of the Coma Cluster, un- 
dertaken using the Advanced Camera for Surveys (ACS). 
We describe the survey and the spectroscopic follow-up in 
section 2 (see also Carter et al. 2008, for further details.). In 
section 3 we present details of our search methods and ex- 
amine a total of seven compact systems which have redshifts 
proving them to be cluster members, three of which are ex- 
treme enough to be classed as cEs. Section 4 provides details 
of our photometric and structural analysis of these objects 
and section 5 an analysis of their spectroscopic properties, 
including velocity dispersions and stellar populations. Sec- 
tion 6 summarises our results. A separate paper will explore 
the range from UCDs to large globular clusters in Coma. 

Before continuing it is prudent to comment on the 
nomenclature used throughout this paper. The term cE is 
used here in a general sense to refer to those galaxies in our 
sample which possess comparable high surface brightness to 
that of M32. It is not used to make assumptions regard- 
ing their specific morphology or indeed possible progenitor 
morphology. 



2 OBSERVATIONS 

2.1 The HST/ ACS Coma Cluster Treasury 
Survey 

The HST Treasury Survey of the Coma Cluster (HST pro- 
gram GO 10861, P.I. D. Carter) currently covers an area of 
just over 200 square arcminutes; unfortunately it was cur- 
tailed when only 28% complete due to the failure of ACS 
on 2007 January 27. The majority of the 25 fields surveyed 
(21 completely, 4 partially) lie in the central regions of the 
cluster, but a few sparsely sample the outer (infall) regions, 
around 1 degree from the core and are centred on known 
cluster members (see ICarter et al.l I2OO8I , for details of the 
pointings). Note that all of the fields observed contain one 
or more moderately bright galaxies from the catalogue of 
(jCodwin. Metcalfe Peachlfl98a CMP). 

Each individual ACS field is 11.3 arcmin^ i n area, im- 
aged onto two 4096 x 2048 pixel CCDs (see iBoffi et al.1 
I2OO7I . for details of the instrument). We assume hereafter 
a luminosity distance to Coma of 100 Mpc (eq uivalent to 
Ho ~ 71 km s~^, see e.g. iKavelaars et al.l [20001 ), i.e. a dis- 
tance modulus m — M = 35.0. At this distance the pixel size 
of the ACS (0.05'') corresponds to 23 pc and the ACS field 
of view is approximately 100 kpc square. Exposure times 
were 4 x 640s in the F475W filter and 4 x 350s in F814W, 
giving limiting point source magnitudes of gF475W = 26.7 
mag, If814w — 25.9 mag (lOcr detections) in the AB sys- 
tem. Standard data reduction was carried out at STScI, us- 



ing PyRAF/STSDAS and MultiDrizzle (|Koekemoer et al.1 
I2OO2I ). Additional se cond-pass pr ocessing is described in the 
main survey paper (ICarter et al.l l2008V Final data products 
and object catalogues are discussed fully in Hammer et al. 
(2009, in preparation), hereafter Paper 2 of this series. 



2.2 MMT / Hectospec Observations 

As part of the ACS Coma Cluster Survey, we are carry- 
ing out a comprehensive spectroscopic survey of galaxies 
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in the Coma refflon usin^ Hectospec, a 3QQ-fibre spectro- 
graph at the MMT (iFabricant et al.ll2005l ). Observing time 
was allocated for this survey by the NOAO TAG through 
the NSF/TSIP program. Hectospec observations were ob- 
tained during three semesters (2007A - 2009A) in a classical- 
queue mode overseen by Nelson Caldwell at the Harvard- 
Smithsonian Center for Astrophysics. We describe the sur- 
vey in detail in Marzke et al. (2009, in preparation). 

Here, we use observations obtained during 2007 and 
2008. For these observing runs, we selected targets drawn 
from SDSS imaging based on their apparent magnitude, 
g — r colour and surface brightness. In the regions observed 
with HST/ACS, the overall spectroscopic completeness is 
90% at r = 19.9 mag and 50% at r = 20.8 mag. Because 
the ACS imaging was unavailable when the first round of 
Hectospec targeting was completed in 2007, we chose to ex- 
clude sources identified as stars by the SDSS photometric 
pipeline. In 2008, we took a similar approach, but we added 
four candidate cEs identified using the procedure described 
in section 3.1. 

We used the 270 lines/mm grating, which delivers ap- 
proximately 4.5 A resolution at a pixel scale of 1.21 A/pixel. 
The useful spectral range is approximately 3700-9000 A. For 
the vast majority of targets (including those discussed in 
this paper) , we integrated for a total of one hour divided into 
three twenty- minute exposures. Seeing ranged from 0.7^^ to 
1.8'^ comparable to the 1.5^^ fibre diameter, which subtends 
roughly 730 pc at the Coma distance. 

We reduced the data using HSRED, an IDL package 
written by Richard Cool and based on the IDL pipeline de- 
veloped for the SDSS. HSRED establishes an approximate 
flux calibration using spectra of bright F stars, which have 
well-measured ugriz photometry and are observed in each 
fibre configuration. 

The HSRED pipeline produces redshifts computed us- 
ing a template- fitting algorithm. We also measured red- 
shifts independently by cross-co rrelating the spectra with 
SDSS templates using XCSAO (jKurtz Minlj Il998l ). We 
selected final redshifts and judged their quality by examin- 
ing each of the 6,259 spectra by eye. 



3 THE SAMPLE 
3.1 Selection Criteria 

To begin our search we first made use of the SExtractor 
(iBertin Arnoutsl[l996l . SE) derived catalogues detailed in 
Paper 2. Based on all the ACS Coma images in both F475W 
and F814W filters this database provided total magnitudes, 
matched aperture colours and first pass size and shape pa- 
rameters for every object in our frames. We next com- 
bined the Hectospec redshift survey results with the imaging 
data and constructed the relevant total magnitude, surface 
brightness and size diagnostic plots to differentiate objects 
of interest that have confirmed Coma membership. To refine 
the search further we impose additional constraints in terms 
of object colour, size and the level of resolution of the image 
('class_star'). 

To constrain object colour we convert the colours of five 
of the known cEs to our ACS filter set by de-reddening and 
K- correcting to rest frame colours where necessary, trans- 
forming to the survey's F475W-F814W colour and then 
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Figure 1. Diagram used to select our initial sample and followup 
spectroscopic targets. Unfilled circles are confirmed cluster mem- 
bers inclusive of the Hectospec 2007 observations. Crosses are 
confirmed background objects. Filled circles show those objects 
that meet our size, colour and star classifier constraints. The data 
presented here is as observed with no corrections for extinction, 
K-correction or cosmological dimming. 



applying Coma reddening and K-correction. Following this 
procedure we obtain a colour cut of 1.1 < F475W-F814W < 
1.5 that can be applied directly to our observed SE photom- 
etry. The cEs used were M32 (V-Ic from HyperLedcQ), NGC 
4486B and NGC 584 6A from the SPSS a nd the two objects 
in Abeh 1689 from iMieske et all (|2005h . Colour transfor- 
mations and K-corrections are derived using the synphot 
routine in IRAF/STSDAS and spectral energy distributions 
(SEDs) from BaSTfl. For the colour transformations we use 
SEDs with -1.27 ^ [Fe/H] ^ 0.4 and 3 Gyr ^ Age ^ 14 
Gyr and for the K-correction we employ a 10 Gyr SED 
with [Fe/H] = -0.66. Tests indicate that the K-correction 
is dependent on metallicity and age of the assumed SED as 
one would expect, with variations of ±20%. This only be- 
comes significant when removing the effects of K-correction 
on the g-i colours of the two cEs in Abell 1689 where it is 
^ ±0.1 mag, and even then it is not sufficient to move either 
galaxy out of the colour cut. Indeed we also find no further 
likely compact candidates within ±0.1 mag of our selected 
colour ran ge. Ga lactic extinction is corrected for according 
to Schlege l et al.1 (|T998). 

Making use of the flux radius parameter provided by 
SE, we impose an upper limit of 500 pc on the radius 
that contains 50% of the flux inside the MAG_AUTO Kron- 
hke aperture. This value was determined by inspecting the 
magnitude-size relation for our data set and largely acts to 
guide the eye to areas of interest in the parameter space. 

Finally we enforce a maximum value for the SE star 
classifier of 0.5 in order to exclude unresolved objects from 
our search. In real terms this cut removes objects with SE 
50% flux radius ^1.5 pixels 35 pc at Coma) and so 
excludes stars, GCs and some UCDs. 

Fig. [T] shows the m— <ij.>r^q diagram used to select 

^ http://leda.univ-lyonl.fr/ 

^ http://albione.oa-teramo.inaf.it/ 
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Table 1. Coma compact galaxy (CcG) photometry and adopted 
redshifts. The full identifier refe rs to the Su rvey visit number in 
which the galaxy is located (see I Carter et al.l |2008) Table 2 for 
details). An additional character, in alphabetical order brightest 
to faintest, is used to differentiate galaxies when a visit contains 
more than one object of interest. 



ID 


Ra2ooo 


Dec2ooo 


ms 


B-Ic 


z 


CcGVl 


13:00:47.67 


+28:05:33.9 


20.72 


2.19 


0.0226 


CcGV9at 


13:00:18.85 


+28:00:33.4 


18.34 


2.04 


0.0207 


CcGV9b* 


13:00:27.31 


+28:00:33.3 


20.71 


2.14 


0.0213 


CcGV12* 


12:59:42.29 


+28:00:55.1 


21.98 


2.17 


0.0261 


CcGVlS 


12:59:59.88 


+27:59:21.6 


21.50 


2.13 


0.0218 


CcGV19a* 


12:59:37.18 


+27:58:19.8 


19.79 


2.19 


0.0254 


CcGVlQb* 


12:59:39.20 


+27:59:54.4 


21.33 


2.07 


0.0236 



the colour cut was extended 0.1 mag blue or redward. Of 
the five probable members four were on the footprint for 
the Hectospec followup run in May 2008 and subsequently 
all four were confirmed as cluster members. Unfortunately 
the fifth could not be observed due to fibre deployment is- 
sues. 

Table [T] lists our sample with SE magnitudes and 
matched MAG_AUTO aperture colours corrected for extinc- 
tion, K-correction and converted to B and Ic bandpasses on 
the Vega system via, 

B = F475WAS + 0.329(F475W-F814W)as + 0.097 
(B - le) = 1.287(F475W-F814W)as + 0.538 

which are derived using the methods described previously. 
The SE Kron-like magnitudes are near total being measured 
within an aperture that contains greater than 90% of the 
light for typical galaxy profiles. Matched aperture photom- 
etry is achieved by running SE in dual image mode. 

In Fig. [2] we present the colour-magnitude relation for 
all confirmed Coma members imaged by the survey. Of our 
compact sample, six of the galaxies sit well above the se- 
quence defined by normal cluster galaxies at their respective 
luminosities with only CcGV9a being close to the general 
trend. It is noted that there is some scatter in the corre- 
lation but the systematic deviation of our sample toward 
redder colours by 0.2 to 0.4 mag is clearly apparent. Put 
another way, almost all the CcGs are ~ 3 magnitudes, or 20 
times, fainter than other cluster galaxies of the same colour. 



3.2 Location in the Cluster 

In Fig. [3] we present a CFHT g band mosaic of the inner 
regions of the Coma cluster overlayed with the positions of 
our sample and relevant giant galaxies for reference. We note 
that throughout this subsection the redshifts and their re- 
spective errors for the bright neighbouring galaxies are taken 
from the SDSS as they were not observed by Hectospec. 

Both CcGVl and CcGV19a are seen in projection close 
to significantly larger galaxies. CcGVl is exceptionally close 
to IC 4045 at a projected distance of ^6.7 kpc and given 
a difference in their radial velocities of /\vr — 117dz53 km 
s~^ it can be assumed they are gravitationally interacting. 
CcGVl 9a is less immediately located near its neighbour at 
^-^23 kpc from NGC 4874, yet still sits in the cDs diffuse 
halo. In this case Avr = 438di65 km s~^. 

CcGV19b and CcGV12 are located -68 kpc and 
101 kpc from NGC 4874 respectively. With the former hav- 
ing At'r = 99±65 km s~^, it is likely to be associated with 
the cD while CcGV12, having Av^ = 650db66 km s"^ with 
respect to it, would appear not to be. 

CcGVlS also has a prominent near neighbour in the 
form of NGC 4886 which is located at a projected distance 
of ^-^28 kpc. Comparing their radial velocities we find /\vr — 
138di54 kms"^ which again is possible evidence for their in- 
teraction. Additionally we note that NGC 4889, the bright- 
est galaxy in the cluster, is a projected distance of —54 kpc 
away with /S.Vr — 35±82 km s""*^. 

Finally, both CcGV9a and CcGV9b are relatively iso- 
lated spatially. CcGV9a has a projected distance of —84 kpc 



"^Spectroscopic membership confirmed following photometric 
selection. tCMP 2777 
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Figure 2. Colour-magnitude relation for all confirmed Coma 
members imaged by our ACS observations. Compact galaxies are 
denoted by triangles. Horizontal lines show the upper and lower 
extent of the colour cut used when selecting the compact galaxy 
sample. 



our initial sample and candidates. Coma membership as- 
signment was based on data from NEE|f| and our 2007 spec- 
troscopic data. The mean surface brightness used here is 
computed as, 

<M>H5o= rnF475W + 2.51og[27ri?io] 

where mpAibw is the SE Kron-like magnitude and R^q the 
radius of a circular aperture which contains 50% of the total 
flux within the Kron-like aperture in arcseconds. 

This process revealed three confirmed Coma compact 
galaxies (CcG) with a mean surface brightness well above 
the global Coma trend (CcGVl, 9a and 18 in Fig. 1). In 
addition it highlighted a number of other candidate objects 
without redshifts which were in turn eyeballed and a quali- 
tative probability of cluster membership assigned. Five were 
deemed to have a high probability of membership based on 
their near circular appearance and lack of visible substruc- 
ture. No new high probability candidates were found when 



http:/ /nedwww. ipac.caltech.edu/ 
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Figure 3. Portion from a CFHT MegaCam g-band image of the central part of the Coma 
cluster. The locations of our compact galaxies are indicated, as well as the three giant galaxies 
with which some of them may be associated. The darkened regions show areas outside of the 
HST/ACS survey footprint. 



to NGC 4889 and few major galaxies within this radius. Sim- 
ilarly CcGV9b has no major galaxies within 55 kpc. 

4 PHOTOMETRIC PROPERTIES 
4.1 Structural Parameters 

We have analysed our sa mple ^sin^ the two-dimensional 
fitting algorithm Galfit (Peng et al. 2002) to recover their 
structural parameters. Each galaxy's light profile is mod- 
elled via axially symmetric ellipses and fit with one or more 
user specified analytic functions. The free parameters of each 
function are adjusted until the residual between the real 
and model image (summed over all pixels) is minimised. As 
well as the imaging data Galfit has a further five important 
inputs which are discussed in detail below. These inputs are 
the uncertainty map, the instrumental point spread function 
(PSF), image masks, the choice of parametric light profile 
model and finally the sky level. 



4.1.1 Uncertainty Maps 

To compute a meaningful for a given model and set of pa- 
rameters Galfit requires a per pixel uncertainty or sigma im- 
age. Such maps may be created by the program itself or oth- 
erwise supplied by the user. Here we opt to take advantage 
of the 'ERR' inverse variance map produced by Multidriz- 
zle which includes poisson noise from both sky and object 
together with instrument dependent contributions due to 
dark current and read noise. The maps are scaled by ^0.77 
to correct for correlated pixel noise (D. Hammer 2008, pri- 
vate communcation) and converted to la per pixel images 
ready for input into Galfit. 



4.1.2 PSF 

A well constructed PSF is of primary importance when fit- 
ting models to the inner regions of galaxy surface brightness 
profiles. Galfit convolves the input PSF with the model cur- 
rently being tested before comparing it to the real image and 
determining a goodness of fit for that model. If the object 
being fit is approaching the resolution limit of the instru- 
ment, as is often the case when attempting to decompose 
our compact galaxies into two components, the dependence 
of the derived structural parameters on the PSF increases. 

It is well known that the PSF varies with position across 
the two ACS/WFC chips and so a separate PSF must be 
created for every desired location on the detector. The stan- 
dard recourse is to either use artificial PSFs generated by 
the Tiny Tim software packagjf] or else some sort of empiri- 
cal approximation. Here we explore both methods to obtain 
suitable input PSFs. 

To create our artificial PSFs at the required position for 
each galaxy we make use of DrizzlyTinfl. Briefly, the task 
transforms the x-y coordinates of each galaxy's centroid on 
the drizzled science frame back to the system of the indi- 
vidual geometrically distorted FLT images. It then executes 
TinyTim with the specified filter and a PSF diameter of 
and injects the distorted output PSFs into blank FLT files. 
Next these files are passed through Multidrizzle using the 
same configuration as used for the science frames to pro- 
duce an undistorted PSF at the observed location of each 



^ See http:// www . st sci . edu / software / 1 iny t im^ 
^ Written by Luc Simard. 
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Figure 4. Compact galaxy images in the F475W filter (left column), residuals from subtracting the 
best fitting seeing convolved single Sersic model (middle column) and overall best fitting model (right 
column). All images are lO^'x 10'' (4.6 kpc x 4.6 kpc). 
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Figure 5. Surface brightness profiles for the sample galaxies. The observed profile and best fitting 
model are plotted as filled circles and dashed lines respectively. For galaxies where the best fitting model 
consists of two Sersic profiles, the inner and outer components are denoted by dotted and dash-dotted 
lines respectively. 
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object. Finally the PSF is cutout and centred according to 
the requirements of Galfit|fl. 

In the derivation of our empirical PSFs it would be 
ideal to use stars observed at the same time as our galaxies. 
Unfortunately owing to the high galactic latitude of Coma 
there are an insufficient number to correctly sample the de- 
tector's PSF variations. As a work around we select recent 
observations of Galactic globular cluster 47 Tucanae taken 
during program 10737 (P.I. J. Mack). The images are re- 
duced using Multidrizzle with the same setup as our sci- 
ence data. The IRAF implementation of Daophot II (|Stetsonl 
Il987l : IStetso~ et al ]|T990) is then used to construct a semi- 
empirical PSF that varies quadratically across the CCDs 
from ^500 stars. The PSF has a Moffat function as its an- 
alytic base with six empirical look-up tables used to correct 
the function to the observed data. To sample the extended 
wings of the ACS/WFC PSF we include non-saturated pix- 
els from centrally saturated stars. After the initial run of 
the task we use the current version of the Daophot PSF to 
subtract all stars used to derive the PSF from the 47 Tuc 
science image. Stars which have a high reduced deter- 
mined by Daophot during the subtraction, and so are likely 
having a detrimental effect during the initial phase of PSF 
creation, are removed from the input star list and the PSF 
derivation repeated. Following the completion of this proce- 
dure we generate a PSF at the precise pixel position of each 
galaxy, again centred correctly for input to Galfit. 

With both PSFs in hand we test their consistency when 
used to derive structural parameters. Results are compara- 
ble when fitting the sample with a one component model. 
However, when fitting a two component model we find sig- 
nificant disagreement in n and Re (see equation (1) below) 
for the majority of our galaxies. To get a handle on which 
derived PSF best represents the actual instrumental PSF 
at the time of observation, we select high S/N stars within 
^4Q" of each target galaxy as new input PSFs. Running 
Galfit once more indicates a broad consistency between pa- 
rameters determined using the empirical and stellar PSFs. 
Based on this result we opt to use the Daophot derived PSFs 
when obtaining structural parameters for our sample. 



4.1.3 Image Masks 

When fitting galaxies which have a near neighbour whose 
Ught distribution does not visibly interfere with the target 
galaxy, but may still have some influence on its structural 
parameters due to its presence in the fitting box, we choose 
to mask the secondary. The masks used are based on the re- 
sults returned by our SExtractor analysis of the Coma imag- 
ing. For each neighbour galaxy we take SExtractor's total 
aperture and scale it linearly by a fac tor of 1.5 (2.25 i n area ) 
in a similar approach to that used bvl Haussler et al. 1(|2003). 
In certain situations where there is a significant blending of 
a target and large nearby galaxy, such as CcGVlQa which 
sits in the halo of NGC 4874, we subtract off the neighbour 
using the IRAF/ELLIPSE task. 



^ Input PSFs for Galfit must be centred correctly, N/2 + 1 for an 
even number of pixels per side and N/2 + 0.5 for an odd number. 



4-1 '4 Light Profile Parametrisation 

Rec ent work by ICrahaml (120021) , fChilingarian et al.1 ^00' 



and lChilingarian &: MamonI (|2008 ) has found cEs to require 
a two component fit to their surface brightness profiles. To 
search for similar structure in our galaxies we fit single and 
double ISersid (|l968 ) models which are of the form. 



I{R) = exp 



l/n 



(1) 



where le is the intensity at the half-light or effective radius 
Re. The parameter n defines the shape of the profile with 
n=l being an exponential profile and n=4 the standard R^^^ 
law f de Vaucouleursiri94 8h. The shape parameter is closely 
connected to k by the expression k, ^ 2n — I . 

The decision to fit a double Sersic model was made 
to avoid unnecessary constraints being placed during two- 
component fitting. Indeed as it is currently unclear as to the 
progenitor of cEs, there would appear to be no reason to 
enforce a classic bulge-disk decomposition. 

To test the resolvability of the inner components which 
Galfit detects, we also make use of the task's ability to fit 
a point source + Sersic model. The advantage here is that 
Galfit uses the input PSF, making no further analytic ap- 
proximations, to model the central component. This model 
describes the case of an unresolved nucleus within a Sersic 
halo. The quality of the resultant fit is assessed based on 
the Galfit computed reduced visual inspection of the 

residual map obtained by subtracting the PSF convolved 
model from the data. 



4.1.5 Sky Background Determination 

Making a robust sky background estimation is of key im- 
portance when fitting observed galaxy surface brightness 
profiles with parametric models (|de Jong||l996l V Using the 
Sersic profile as an example, if the sky level is under esti- 
mated then the galaxy being fit appears to have more flux in 
its outer regions. This results in the model's wings becom- 
ing larger and either an artificially inflated n value, a larger 
effective radius or both. If the sky is over estimated the re- 
verse occurs. An additional complication comes from sky 
gradients generated by diffuse intracluster light. It is there- 
fore clear that a reliable means to determine each galaxy's 
local sky background is critical when obtaining their struc- 
tural parameters. 

Galfit allows the sky to be fit as an extra parameter but 
this relies on the assumption that the model profile being 
used well describes the real light profile. Tests indicated that 
this assumption was adequate for one component fits, how- 
ever, when employing a two component model it was found 
that one of the profiles was able to suppress the sky and 
consequently grow too large in both n and Re. To solve this 
problem we determine an d fix the sky pri o r to p rofile fitting 
using a method similar to iHaussler et aD (|2007f ). 

Each galaxy's light distribution is measured by the 
IRAF/ELLIPSE task with a maximum semi-major axis for 
iterative fitting set to the SExtractor total aperture radius. 
In this way the annuli beyond this point are of fixed ellip- 
ticity and orientation. The profile gradient between these 
annuli is computed and the sky value determined as the 
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Table 2. Structural Parameters for CcGs. Col. 1: Galaxy identifier. Col. 2: Best fitting model Sersic + Sersic (S+S) or point source 
+ Sersic (PSF+S). Col. 3-10: Absolute magnitude and mean surface brightness within the effective radius Re in the B band (Vega), 
Sersic n and effective radius Re for inner and outer components respectively. . 



CCG 


Model 


MB,in 


<MS >e,in 




Re, in 


MB,out 


<MS >e,out 


^out 


^e,out 


(ID) 




(mag) 


(mag arcsec"^) 




(pc) 


(mag) 


(mag arcsec"^) 




(pc) 


VI 


PSF+S 


-8.51+0.13 








-14.07+0.01 


21.42+0.01 


1.04+0.01 


252+2 


V9a 


S+S 


-15.72+0.21 


18.84+0.22 


3.49+0.23 


172+36 


-16.22+0.13 


20.95+0.13 


1.64+0.07 


575+10 


V9b 


S+S 


-12.72+0.24 


21.26+0.25 


2.95+0.26 


145+33 


-14.05+0.07 


22.27+0.07 


1.17+0.03 


359+3 


V12 


PSF+S 


-10.42+0.04 








-13.03+0.01 


21.39+0.01 


1.88+0.02 


158+2 


V18 


PSF+S 


-10.04+0.03 








-13.59+0.01 


21.48+0.01 


1.34+0.01 


201+1 


V19a 


S+S 


-14.07+0.09 


18.46+0.1 


4.38+0.19 


76+9 


-14.89+0.04 


20.63+0.04 


1.37+0.03 


254+2 


V19b 


S+S 


-12.90+0.20 


16.82+0.21 


2.60+0.52 


23+5 


-13.19+0.10 


21.09+0.13 


1.40+0.04 


142+14 



mean of the annuli outward of the point at which the gra- 
dient reaches a prescribed small value. The radius at which 
this transition occurs is scaled by a factor 1.5 and results 
in a region of width 25-65 pixels (area ^-^9500-65000 pixels^) 
being used to estimate the sky, depending on galaxy size. 

The above procedure makes no attempt to measure sky 
gradients. Galfit, when permitted to fit the sky freely, mod- 
els it as a plane which may tilt. We take advantage of this 
functionality and allow the task to measure the plane's tilt 
during fitting, while keeping the mean sky level fixed. 

4-1.6 Results 

Table [2] presents the best fitting structural parameters for 
all seven CcGs. Native F475W magnitudes are corrected for 
Galactic extinction and K-correction before being converted 
to the B band using the same methods as before. The mean 
surface brightness parameters are additionally corrected for 
cosmological dimming. All magnitudes and surface bright- 
ness measurements are converted to a Vega zeropoint|j for 
comparison with previous work. 

In Fig.Uwe present three cutout images of each galaxy. 
The left hand column is the observed data, the middle col- 
umn is the residuals produced by subtracting the best fit- 
ting seeing convolved single Sersic model from the data and 
the right column contains the equivalent but for the model 
which best describes each galaxy (see Table [2] column 2) . 
To further display the quality of the best fits in Fig. [5] we 
use IRAF/ELLIPSE to produce surface brightness profiles of 
both real and model data and hence obtain one- dimensional 
residuals. 

It is evident that the majority of our sample are well fit 
by their respective models. All galaxies were found to sta- 
tistically prefer two components, either two Sersic profiles 
(S+S) or a point source plus Sersic profile (PSF+S), over a 
single Sersic model. Only CcGVQa and CcGV12 have resid- 
ual images showing remaining substructure. Based on our 
colour maps constructed in the next section, it is unlikely 
the residuals associated with CcGV12 are part of the galaxy 
given their bluer colour with respect to the rest of the galaxy 
at such radii. CcGVQa residuals are interestingly symmetric 
and are probably indicative of a weak bar. 

Both CcGVl and CcGVlS have comparable reduced 
when fit by PSF+S and S+S models. This likely indi- 
cates the inner component is almost resolved but is either 
too faint, too small or both to be modelled clearly. As such 

^ http:/ /www. stsci.edu/hst/acs/analysis/zeropoints 



we opt for the simpler PSF+S model to parameterise these 
galaxies. Also of note is the inner component of CcGVlQb 
which has an Re = 23 pc or 1 ACS/WFC pixel. Given 
its exceptionally small angular size both resolvability and 
repeatability checks were conducted using Galfit. 

Fitting a PSF+S model and comparing this with a S+S 
fit we find Axl — 0.11 in favour of the latter. We also per- 
formed some simulations using Galfit to first create a set of 
artificial S+S models with 0.2 ^ Re,in ^1.0 pixel in steps 
of 0.05 pixel, convolved with the PSF used for CcGVlQb. 
Next we fitted the models with both PSF+S and S+S pro- 
files and estimate our limit for detection of a resolved inner 
component for this particular galaxy to be Re,in ^ 0.4 pixel 
or 9 pc. From this we conclude that the core of CcGVlQb is 
sufficiently well resolved, likely due to its exceptional lumi- 
nosity. We also note that the inner component of this galaxy 
is equally well fit by a King profile, often used when studying 
GCs. Fitting King + Sersic we find that the derived struc- 
tural parameters for both inner and outer components are 
consistent to well within the quoted errors in Table (2] 

Typically Galfit estimates errors base d on the shape of 
the surface surrounding the Xmin {see Peng et al.l l2002l 
for further details) and in general provides reasonable uncer- 
tainty estimates. However in certain regimes there are likely 
to be further systematics which may become more significant 
such as in the case of CcGVlQb where the observed profile is 
poorly sampled and CcGVl where the point source compo- 
nent is extremely faint. By varying the starting values of Re 
and n given to Galfit we check the stability of our derived 
parameters for the inner component of CcGVl 9b, which are 
reported by Galfit to have uncertainties of around a few per- 
cent. Such tests indicate a more realistic uncertainty of ^ 
20% for both Re and n while the component's magnitude 
varies by roughly ± 0.2 mag. For CcGVl, we use the visit's 
uncertainty map to produce new realisations of the frame 
and then pass them through Galfit. The standard deviation 
of the distribution of magnitudes for the point source com- 
ponent is then used as its quoted error in Table [2] 

To place our sample in context with other recent work in 
this area we present in F ig.[6]updated versi ons of Figs. 9a g 
fromlGraham Guzma n (2003) (see also .Chilingarian et alj 
l2007l ). These diagrams compare mean surface brightness 
within the effective radius Re with absolute magnitude in 
the B band and effective radius Re for dynamically hot stel- 
lar systems. They represent differe nt versions of the well 
known Mb - fJ^ef f iKormendvl (| 19771 ) re lation which itself is 
a projection of the Fundamental Plane (|Diorgovski DavisI 
119871 ). For the full range of elliptical galaxies, dwarf 
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Figure 6. Structural parameters for dwarf to giant elliptical galaxies, bulges, nuclei and compact galaxies from the literature and this 
work (purple filled squares). The left panel displays the parameters for the combined or whole profile of all CcGs, M32 and A496cE while 
the right panel presents the inner components of our sample where they are resolved, M32 and A496cE. Coma compact galaxy points are 
not labelled in some panels for clarity's sake. 



to giant, d a ta is t aken from iBinggeli JerienI (Il998l). 
ICaon et alJ (Il993l\ I P'Onofrio et al.' ('1994V 'Faber et al.' 
(|l997l l. iGraham k G uzman (2003), Stiavelli et al. (2001) 
and homogenised as per Grah am Guzmanl ()2003h . We ad- 
ditionally inc lude data for Es and d Es from the Virgo Cluster 
ACS Survey (|Ferrarese et al.|[2006h , with the gF475W photo- 
metric parameters transformed to the B band using sim- 
ilar methods to those described for ou r Coma data, and 
Es brighter than Mb = -18 mag from |Prug;niel Simien 
( 1996i) . A sample of bulges is taken from iMac Arthur et al. 
( 20031 ) . In ter ms of compact galaxies from t he literature we 
include M32 (Graham" '2002') , NGC 4486B (jFerrarese et al.1 
[2006) and A496cE (Chilingar ian et al.ll2007h . 

It is current practice to assume cEs that are best fit by 
two component models such as Sersic + Exponential disk or 
S+S with outer n ~ 1 are indeed bulge/disk decompositions. 
As such only their bulges would be plotted on diagrams like 
Fig.[6l That said, quite whether these outer components are 
in fact rotationally supported disks or just the remains of 
a dwarf/intermediate- luminosity elliptical progenitor galaxy 
are unclear. With this in mind we present two versions of 
both plots. The left panel of Fig. [6] shows the position in 
the parameter space of the combined two component mod- 
els for ah CcGs, M32 and A496cE. In this case the "whole" 
profile Re is obtained via numerical integration. The right 
panel of Fig. [6] then presents the structural parameters for 
the Sersic inner components of all Coma compacts where 
these are resolved, M32 and A496cE. The position of NGC 
4486B does not change in these plots as it is modelled by 
a single Sersic function in the literature. To aid comparison 
we also plot data o btaine d for Fornax and Virgo UCDs by 
lEvstigneeva et aL (l2007b) and galaxy nuclei in the Virgo 
Cluster from ICote et al. 1 120061 ). For those UCDs which are 
best fit by an inner compact core and outer halo, we plot only 
the inner component. Here we note that both the UCDs and 



the galaxy nuclei have had their structural parameters de- 
rived from King profile fits rather than a Sersic model. From 
the tests conducted for the inner component of CcGV19b 
a nd by inspection of the resu lts presented in Tables 8 and 9 
of lEvstigneeva et"al] (|2007bl ), we typically expect at most a 
~ 30% difference in parameters derived from the two fitting 
functions (VUCD 3 excluded). 

The combined profile diagram clearly shows that three 
of the CcGs, V9a, V19a and V19b, have structural param- 
eters which are comparable, if marginally less extreme, to 
those of confirmed cEs. They are seen to extend the sequence 
defined by giant ellipticals toward fainter magnitudes/higher 
surface brightness in line with those compacts from the lit- 
erature. The remaining four CcGs appear to follow a com- 
parable sequence, albeit with lower surface brightness, orig- 
inating from the dE region of the diagram. 

Turning to the right panel, the inner components of 
CcGV9a and CcGV19a are seen to have similar structural 
properties to the inner components of A496cE and M32 re- 
spectively. The most extreme case is CcGV19b whose posi- 
tion on the diagram is well within the parameter space oc- 
cupied by the more luminous/more massive nuclei of early- 
type galaxies in Virgo. CcGV9b on the other hand appears 
to have an inner component much less extreme in terms of 
surface brightness and relies instead on small physical size 
and faint magnitude to perturb it off the correlation seen for 
normal dEs. It is evident that the inner components of the 
Coma cEs, in combination with those from the literature, 
begin to fill up the region of the plot between Es/bulges and 
UCDs/nuclei. Additionally, given the observed sequence it 
is possible to speculate that, structurally speaking, the inner 
component of CcGV19b, and perhaps CcGV19a, are likely 
galaxy nuclei while that of CcGV9a appears more bulge like. 
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Figure 7. Colourmaps in F475W-F814W filters from the ACS imaging. We also plot each galaxies colour 
and ellipticity profile. 



4.2 Colour Maps 

Colour gradients in early-type galaxies can be direct indi- 
cators as to their formation and evolution mechanisms. If 
one neglects the influence of any dust present in a given 
galaxy, such gradients are generated by variations in stellar 
population age and metallicity with radius. In terms of for- 



mation, monolithic collapse has been predicted to es tablish 
a steep negative metallicity gradient (Carlberd lioii l which 
is then progressively diminished to varying degre es by major 
and minor mergers during a galaxy's evolution ([Kobavashil 
1 20041 ) . Hence colour gradients in early- type galaxies are typ- 
ically considered to be due to metallicity variations alone 
rather than age although the validity of this assumption is 
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still s omewhat in question (jSilva Elstonlll994l :l Baes et alJ 
I2OO7I ). 

Previously M32 has been found to hav e a flat colour 
profile in both the optical and near-infrared (|Peletierlll993l ) 
which is, apparently, du e to an exact com pensation in metal- 
licity and age gradients ("Rose et al ]l2005|).Both A496cE and 
NGC 4486B have also been found to have no colour gradi- 
ent, although given the scale used t o display the colour axis 
of Fig. 81 in IPerrarese eraP (|2006h it is diflicult to assess 
the possibility for a moderate gradient as is often the case 
in early- type galaxies . On t he other hand M59cO discovered 
by Chilin garian & Mamo was observed to have a 

relatively strong gradient in the sense that the galaxy's core 
was ~0.15 mag bluer than its halo. 

To search for radial stellar population gradients in 
our sample we have constructed colour maps and pro- 
files from our ACS imaging. In order to do this we must 
first degrade the resoluti on of the P475W imape to that 
of the P814W. Pollowing Chilingarian & Mamon (2008) we 
smooth the P475W by the ratio of the fourier transforms of 
PSPf814w and PSPf475 using the IRAP task PSPMATCH. 
This method acts to minimise resolution loss which would 
occur when convolving either science image with the PSP 
of the other. We then run IRAP/ELLIPSE on both images, 
taking the isophote table generated from the P475W run 
and using them as input isophotes for the P814W image. 
This process correctly ensures that our colour profiles are 
derived from matched annuli in both filters. In Pig. [71 we 
present the colour profiles and maps constructed using this 
method as well as isophotal ellipticity profiles obtained from 
the P475W images. 

Nearly fiat profiles are observed for CcGVl and 
CcGV19b which may point to similar stellar populations 
at all radii although this is not definitive given the age- 
metallicity degeneracy suffered by broadband colours (and 
the example of M32 noted above). CcGV18 has a shallow 
gradient, the colour becoming ^ 0.2 mag redder as one 
moves inward, though interestingly within the central 0.2^^ 
there is a noticeable blueing in colour of order ~ 0.1 mag, 
perhaps indicating a younger population present in its unre- 
solved nucleus (as in dEs such as NGC205). The remaining 
four galaxies, which include the cEs CcGV9a and CcGV19a, 
all contain strong colour gradients with bluer colours at 
larger radii. 

The colour map of CcGV9a clearly shows a red core re- 
gion with higher ellipticity isophotes than those outward of 
2" . The peak in the ellipticity profile at ^ 1.^" provides fur- 
ther evidence for the existence of a weak bar in this galaxy. 
There is a noticeable change in the colour profile slope at 
^ 0.3^^ which corresponds well with the effective radius of 
the inner Sersic component returned by Galfit. Beyond ^ 
2.2^^ the profile turns over and gradually becomes redder by 
0.15 mag. Prom the colour map alone this galaxy certainly 
appears to be at least a two component system with a signif- 
icant colour difference between the inner and outer regions. 

CcGV19a is observed to have a colour variation of ^ 
0.3 mag within just 2" . As with CcGV9a there is again a 
change, albeit poorly sampled, in colour profile slope near 
the effective radius of the inner component detected by Gal- 
fit, ~ 0.1 to 0.1 Beyond this transition the colour profile 
decreases in a remarkably linear fashion until ^ 1.8^^ while 



the isophotes maintain a fairly stable ellipticity e = 0.1 ± 
0.05. 

The colour profile of CcGV9b has a slope change at a 
smaller radius than would be suggested by its inner compo- 
nent effective radius. This is perhaps due to the smaller dif- 
ference between inner and outer component surface bright- 
ness, A/i = 1 mag for CcGV9b but, Afi>2 mag for CcGV9a 
and 19a, where A/i is the difference between inner and outer 
component </is >e- This aside, CcGV9b and 19a have very 
similar colour gradients, both becoming ^0.3 mag bluer 
outwards, which may be evidence for comparable changes 
in stellar population with radius between the two. Whilst 
speculative, if this were true it would be an interesting re- 
sult given the differences in their structural parameters. 

Pinally CcGV12 also contains a strong colour gradient 
which is mostly confined to within 0.2" . The structure iden- 
tified in the Galfit residual is clearly visible in the colour 
map and is masked when the profile is constructed, using a 
simple binary mask of approximately comparable size. Once 
outside the nuclear regions the colour profile steadily moves 
to bluer colours with increasing radius. As with all the galax- 
ies with strong gradients, it is likely this galaxy contains a 
central component which is older or more metal rich, or in- 
deed both, relative to the outer halo. 



5 SPECTROSCOPIC ANALYSIS 

5.1 Velocity Dispersions and Masses 

To measure the velocity dispersions of our sample we em- 
pl oy the Penaliz ed Pixel Pitting (PPXP) routine developed 
bv fcaDDellaTrfc EmsellemI (|2004l ). This technique works by 
fitting a set of template spectra convolved with a paramet- 
ric line of sight velocity dispersion (LOSVD) to the observed 
galaxy spectrum in the pixel space. Both the target galaxy 
and template spectra are first rebinned to a linear scale in 
In (wavelength) and then a model M of the target's spectrum 
G is constructed such that, 

k I 
M^PmY,^^ (^(^' V^) * Tfc) + ^ &^P^ (2) 

where Pm and Pi are multiplicative and additive Legendre 
polynomials respectively of order m and / which are used to 
correct the continuum shape of the templates during the fit, 
Tfc is a set of template spectra, B{a^V) is the broadening 
function used to parametrise the LOSVD and Wk and bi 
the weights given to each template. Pinding the best fitting 
parameters for B is then a matter of minimising, 

n = l ^ ^ 

where AC is the per pixel uncertainty on G obtained during 
reduction. W hile it is possible to fit for hi g her moments o f 
the LOSVD (jvan der Marel Pranxl[l993l : lGerhardlll993l l. 
given the intermediate velocity resolution of our data we 
choose to fit for a alone setting B as a gaussian func- 
tion. Por Tfc we use the simple stellar population (SSP) 
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models of Vazdekis et al. (in preparatior0) which are 
based on the MILES empiri cal stellar spectral library 
{'Sanchez-BlazQuez et al." 2006a). These were chosen due to 
their large spectral coverage of 3540-7410A and intermedi- 
ate resolution of 2.3A FWHM which is better (by a factor of 
2) than that of our data, as required if they are to be used 
as templates in PPXF. To determine the order of Pm and 
we fix / = 2 and increment m — 0-10 in steps of unity. 
We then select m based on a trade off between the stability 
of a and the increasing computational time required as the 
complexity of Pm grows. Typical values are m = 3-6 and for 
Pi we find that / = 2 satisfies our criteria. 

The actual fitting proceeds by smoothing the template 
spectra to match the instrumental resolution and rebinning 
data, templates and uncertainty spectra onto a logarithmic 
scale. Next the fitting range is set to 4000-6500A (in the 
rest frame), a region chosen so as to avoid the saturated H 
and K lines in the blue and the increasing sky subtraction 
uncertainties in the red. Taking advantage of fitting in pixel 
as opposed to Fourier space, we mask poorly subtracted sky 
lines and the NaD resonance line at 5892 A. As our spec- 
tra were taken only with measuring redshifts in mind, their 
signal- to- noise (S/N) is often less than ideal. However, we 
make an attempt to fit the entire sample using PPXF and 
to suitably quantify the uncertainty on each LOSVD for 
completeness. Fig. [8] and Table [3] present the spectra of our 
sample together with best fitting model, residuals of data 
minus model and la uncertainty spectra and the velocity 
dispersion resulting from the fits. 

We test the robustness of our measurements in a variety 
of ways. To estimate the random error present on our veloc- 
ity dispersions we perturb the best fitting model spectrum 
for each galaxy by the relevant uncertainty spectrum and 
run PPXF again on the new pseudo-observation. Following 
100 executions of this process we take the standard deviation 
of the distribution of velocity dispersions as our la errors 
quoted in Table O As we use the best fitting model for this 
process, template mismatch errors are not propagated into 
this uncertainty. 

Systematics in terms of the routine's ability to resolve 
dispersions significantly below our velocity resolution, par- 
ticularly from data of low S/N, are again tested for via sim- 
ulations. We take a high S/N model spectrum of comparable 
metallicity and age as our galaxies (see next section) from 
the Vazdekis library and smooth it to our instrumental res- 
olution. Next we convolve the spectrum with a gaussian of 
width equal to the measured velocity dispersion of a particu- 
lar galaxy in our sample and add the relevant noise. Starting 
with CcGV12, CcGVlS and CcGV9b we repeatedly gener- 
ate such fake observations with known velocity dispersion 
and run them through PPXF. These galaxies possess the 
lowest S/N and smallest measured velocity dispersions of 
our sample. We find CcGV12 and CcGVlS to have system- 
atic uncertainties in their velocity dispersions of 10% while 
CcGVQb has 25% uncertainty on its measurement. To place 
this in context, CcGVl has a ^1% systematic error. 

As our wavelength fitting range is larger than typi- 
cally used when measuring velocity dispersions, we also test 
PPXF over a range of 4400-5400A. In ah cases we find that 
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Figure 9. Faber-Jackson relation for elliptical galaxies, UCDs, 
DGTOs, GCs and literature compact galaxies. Coma compacts 
are plotted as black filled circles. The horizontal arrow ends show 
the absolute luminosity of the central components of our sample 
where resolved. The vertical arrow ends correspond to the central 
velocity dispersion measurements obtaine d for M32, NGC 4486B 
and NGC 5846A bv lDavidge et aT 



the differences are consistent within the random error pre- 
sented in Table [3] Therefore, along with the resolvability is- 
sues discussed above, we make no attempt to correct our ve- 
locity dispersions for systematics as it is clearly the random 
error that dictates the uncertainty in our measurements. Fi- 
nally we note that CcGV9a is the only galaxy in our sam- 
ple that has had its velocity dispersion measured elsewhere. 
SDSS DR6 observed the galaxy as having a velocity disper- 
sion of 87±9 km s~^ through a 3^^ diameter fibre which is 
consistent with our result wit hin the joint unce r tainti es. 

In Fig. [9] we present the iFaber JacksonI (|l976l l rela- 
tion for a compilation of h ot stellar systems. For elliptical 
galaxies w e take data from Faber et al. (1989), Geha et al.J 
(200 3). van Zee et all 12004) Matkovic k Guzman (2005) 
and IChihngarian et al.n2008l). We p lot data for GCs in 
NGC 5128 from iMartini Hoi (12 004*) and the Milky Way 
from iMcLaughlin van der Mar el (2005) as w ell as UCDs 
i n For nax from Drinkwater et al. (2003) and Hil ker et al.l 
(|2007l ) and DGTOs in Virgo from Hasegan et al. (20051 ). Fi- 
nally we include se veral cEs from the l iterature. For M32 
we t ake data from 'Be nder et aP (|l996l ). for NGC 4486B 
fr om|Smith eT al. (2000), for NGC 5846A from|Muller et al 



1999) and the SDSS, fo r A496cE from IChihngarian et al 



(2007) and M59cO from IChilingarian &; Mamonl (|2008l ). B 
band absolute magnitudes f or M32 and NGC 4486B are 
those reported in Table 1 of 'Chilingarian et al.1 (|2007l ) and 
references therein. It should be noted that the three archety- 
pal compact galaxies, M32, NGC 4486B and NGC 5846A, 
have had their central velocity dispersions measured a num- 
ber of times in the literature with a variety of differ- 
ent instruments. Using HyperLeda we have selected up to 
date representative measurements from ground based stud- 
ies which provide a suitable comparison to the data pre- 
sented in the rest of the diagram. Re cent high spatial reso - 
lution adaptive optics observations bv lDavidge et all (|2008l ) 
of these three galaxies are also plotted for completeness and 
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Figure 8. CcG spectra as observed by Hectospec. Observed data, PPXF best fits, residuals and la error spectra are denoted by blue, 
red, cyan and black lines respectively. 
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Figure 8. continued. 



correspond to the upper ends of the arrows attached t o thes e 
galaxies. Using the prescription of |j0rgensen et aP (|l995h . 
we correct our velocity dispersion measurements to a suit- 
ably small aperture (total i?e/4) to estimate their central 
values. In all cases the correction is small, < 13%, and com- 
parable to the measurement errors. The estimated central 
velocity dispersions for our sample are tabulated in Table [3] 
and plotted on Fig. [9] The head of each arrow denotes the 
position on the diagram of the central component of each cE 
where resolved. For M32 and A496cE only the bulge abso- 
lute magnitude is plotted. We note that conventional galax- 
ies are largely absent from the parameter space of Fig. [9] at 



Mb < -15 mag, likely owing to their low surface brightness 
(see Fig. [6]). 

Fig. [9] shows that, along with the cEs from the litera- 
ture, the majority of compact galaxies presented here con- 
form to a tight sequence beginning at DGTOs and UCDs 
and extending to well above normal elliptical galaxies. Both 
CcGV9a and CcGV19a have velocity dispersions that are 
a factor of two larger than dEs at their respective magni- 
tudes. It is noted that intermediate luminosity ellipticals do 
almost reach CcGV9a/A496cE, though they are clearly ex- 
treme outliers. The Coma compacts, CcGV9b aside, essen- 
tially begin to fill in the parameter space between M59cO 
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Table 3. Median S/N per A and measured velocity dispersions for our sample based on Hectospec 
2007 and 2008 observations and a fitting range of 4000-6500A. The fourth and fifth columns present the 
estimated central and global velocity dispersions for each galaxy while the final two columns display the 
virial mass estimate and corresponding M/L within 5Re. 



CcG 


S/N 


cr 


cro 


(^global 


^dyn 


{M/LB)dyn 


(ID) 




(km s-i) 


(km s-i) 


(km s-i) 


(108 Mo) 


(Mo/Lo) 


VI 


13.2 


85.2±13.9 


91.2±14.9 


80.9±13.2 


37.9±12.4 


57.1±18.6 


V9a 


38.4 


97.3± 5.0 


102.3± 5.4 


90.7± 4.8 


75.6± 7.9 


9.6± 1.4 


V9b 


12.7 


25.1±15.8 


26.6±16.8 


23.6±14.9 


4.0± 5.1 


4.8± 6.1 


V12 


7.3 


42.2±30.0 


46.4±33.0 


41.2±29.3 


5.0± 7.1 


18.1±25.7 


V18 


7.1 


66.3±25.3 


71.9±27.4 


63.8±24.3 


17.3±13.2 


45.3±34.5 


V19a 


33.9 


76.1± 5.1 


82.3± 5.6 


73.0± 4.9 


24.2± 3.3 


11. 7± 1.6 


V19b 


13.9 


65.8±17.7 


74.1±19.9 


65.7±17.6 


7.2± 3.8 


13.8± 7.5 



and M32. From the diagram only CcGVQb appears to differ 
from the trend. Unfortunately due to the large uncertainty 
on this galaxy's velocity dispersion it is hard to comment in 
a definitive way. However it would seem possible that this 
galaxy is some sort of intermediate between dEs and the 
compact sequence. 

With the velocity dis persions in han d we are able to use 
the virial mass estimator (ISDitzerlllQSTh . 



G 



(4) 



where a is the global projected velocity dispersion, Re the 
projected half light (mass) radius and [3 the virial coefficient. 
We estimate a by correcting our data via the expression of 
|j0rgensen et al. (1995) to an aperture of radius 5i?e and 
tabulate our estimated global dispersions in Table [S] For all 
profiles considered this contains >90% of the flux and by 
inference mass of the system. Total galaxy Re is again com- 
puted by numerical integration. The choice of scaling factor, 
here set as /3 = 10, warrants some further comment. This fac- 
tor incorporates a number of corrections including those for 
projection eflFects and the specific d ensity distribution of th e 
system in question. Recent work bv lCapDellari et all (|2006h . 
which compared complex dynamical models with equation 
(4), found that for a sample of 25 E/SO galaxies /3 '-^ 5 with 
cr measured within an aperture of radius Re. This result 
is specifically based on one component R^^^ profiles and is 
demonstrated in their Fig. 13, which shows the relation be- 
tween their computed masses and those from equation (4) 
assuming /3 = 5. It is interesting to note that in this figure 
M32, the lowest mass galaxy in their sample, is a substan- 
tial outlier above the relation, well beyond quoted errors. 
Disregarding other cEs, the next step down in the compact 
stellar system hierarchy are UCDs. Studies of these objects 
in both Fornax and Virgo have found that /3 ^ 10 repro- 
duce s well the results of more ri gorous methods (|Hilker et al.l 
I2OO7I : lEvstigneeva et al.ll2007bh . Therefore, as it would seem 
(3 could in fact be anywhere between these two extremes for 
our obj ects, we opt to follow pr e vious compact galaxy work 
such as IChilingarian MamonI (|2008[ ) and set (3 — 10. 

Dynamical masses and mass to light ratios (M/L) for 
our galaxies can be found in Table [S] Uncertainties are those 
obtained by propagating the error on each galaxy's velocity 
dispersion and total B band magnitude. Aside from random 
errors, there are a number of caveats which apply when using 
the virial mass estimator. The expression assumes a constant 
M/L with radius, an isotropic velocity distribution and that 



the galaxies are spherical. Given these factors and the inher- 
ent uncertainty in scaling the aperture velocity dispersions 
to global values, the masses and corresponding M/L pre- 
sented in Table [3] should be treated as best estimates given 
the available data. 



5.2 Stellar Population Parameters 

The standard approach when attempting to recover lumi- 
nosity weighted age and metallicity for unresolved stellar 
po pulations is that of li ne in dex measurements. Pioneered 
bv lBurstein et al.1 (|l984l ) and lGorgas et all (|l993h these in- 
dices target key information rich absorption features present 
in integrated galaxy spectra . Following the introduction of 
Lick Indices in particular (F aber et al.l 19851 ), it was shown 



that when combined with predictions from SSP models they 
are able to break the age- metallicity degeneracy for such sys- 
tems inherent in broad band photometric models (|Worthevl 
"1994"). More recent work has focused on updating and refin- 
ing SSP models to corr ectly quantif y the effect s of non -sola r 
abundance patterns (| Thomas et al.l [2OO3I : ISchiavonl [200^ 
whic h have been observed in early-type galaxies for some 
time (IPelet iedl 19891 : 1 Wort hev et al.ll992l : lTragei][i997h . Typ- 
ically such a-enhancement has been detected as an over 
abundance of Mg relative to Fe in luminous elliptical galax- 
ies and has been seen as evidence for shorter star- format ion 
episodes in more massive galaxies. This stems from the en- 
richment history of the interstellar medium out of which 
the current generation of stars formed. On short timescales 
this is dominated by Type II supernova, yielding more Mg 
than Fe relative to the solar neighbourhood, but on longer 
timescales by Type la supernova which red ress the balance . 

In this section we use the models of ISchiavon ( 2007 ) 
and the EZ-Ages code devised bv lGrave"s SchiavonI (|2008h 
to measure the age, metallicity and key abundance indi- 
cator [Mg/Fe] of our compact galaxies. Brie fly, th e Schi - 
avon models are based on the flux-calibrated Ijonesi (|l999h 
spectral library and allow for variable abundance ratios in 
terms of Mg, C, N and Ca to be present in the stellar atmo- 
spheres. Additionally the model spans a range of 1 to 17.7 
Gyr and -1.3 to +0.2 in [Fe/H]. Note that the models are 
cast at given [Fe/H], rather than total metallicity Z. Here we 
use a Salpeter initial stellar mass function and scaled-solar 
isochrones. 

EZ-Ages employs the model via a sequential grid inver- 
sion method to return consistent stellar population param- 
eters. The process begins by making an initial measurement 
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Table 4. Absorption line index data used in our analysis of the 
CcGs stellar population parameters. 



CcG 


Hbeta 


Mgb 


Fe5270 


Fe5335 


(ID) 


(A) 


(A) 


(A) 


(A) 


VI 


1.94±0.45 


3.64±0.44 


3.30±0.46 


3.23±0.52 


V9a 


1.74±0.14 


4.09±0.16 


2.16±0.25 


1.86±0.19 


V9b 


2.62±0.45 


3.72±0.45 


2.67±0.48 


2.23±0.56 


V19a 


1.62±0.17 


4.50±0.17 


3.00±0.18 


2.41±0.21 


V19b 


1.71±0.42 


3.88±0.40 


1.87±0.45 


2.09±0.53 



of age and metallicity from the Lick H/3 and <Fe> indices 
of the galaxy in question, where <Fe> is the average of the 
Fe5270 and Fe5335 indices. These are selected as they pro- 
vide good age (H/3) and metallicity response (Fe5335 and 
Fe5270), while bei ng; essentially i nsensitive to abundance 
pattern variations ([Schiavonl l2007l ) . From the HP - <Fe> 
index-index grid a two dimensional linear interpolation con- 
verts the position in index space to the fiducial age and 
metallicity values. 

Next the code moves onto determining the [Mg/Fe] ra- 
tio and to do this creates a new grid of H/3-Mgb5177 (or 
Mg2). If the age and metallicity obtained from the inversion 
of this grid do not match the initial estimates it is assumed 
a non-solar [Mg/Fe] is present and the model is recomputed 
with a new [Mg/Fe] ratio. Further iterations occur until the 
age and metallicity derived from both grids match within 
some tolerance. EZ-Ages then goes on to repeat this pro- 
cess replacing Mg for C, N and Ca in turn with the relevant 
index used to compare back to the fiducial age and metal- 
licity estimates. Finally the new abundance pattern is used 
to compute a new model and derive a new fiducial age and 
metallicity with the entire fitting sequence repeated until 
the fit does not improve. 

Unfortunately the S/N of our spectra is once again an 
issue and we restrict our line index analysis to the five galax- 
ies with the highest quality spectra to avoid uncertainties 
that would span entire index-index diagrams. Furthermore, 
while a full abundance pattern is fit, as the uncertainties as- 
sociated with the initial parameters are propagated in turn 
through the fitting scheme, we limit our analysis to age, 
metalUcity and [Mg/Fe]. 

To measure the desired indices we use the IDL code pro- 
vided with EZ-Ages. In addition to the galaxy and uncer- 
tainty spectra the script requires the velocity dispersion for 
each galaxy so that the spectra may be correctly smoothed 
to the Lick reso lution which itself varies for each index (see 
I Schiavonl (|2007l l Table 1). The required smoothing in each 



case IS (T smooth 



j2 

lick 



res 



gal 



where ares is the 



Hectospec instrumental resolution and agai the galaxy ve- 
locity di spersion. Inde x error s are computed using the equa- 
tions of ICardiel et all (|l998l l. The relevant indices used in 
our analysis are tabulated in Table HI 

Fig. [To] presents an example H/3-<Fe> diagram used 
by EZ-Ages to determine fiducial ages and metallicities. It 
is important to note that the index grid in this figure is gen- 
erated from a model with a solar abundance pattern and 
so is slightly offset with respect to the grids that would be 
obtained from each galaxy's final best fitting abundance pat- 
tern. Hence to avoid overlaying five only marginally different 
grids we opt to display the data this way for reference only. 

Table [5] displays the stellar population parameters of 
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Figure 10. Example of the index-index diagram used by EZ- 
Ages to measure fiducial age and metallicity before going on to 
fit a full abundance pattern. The grid is for a solar abundance 
model and so is not tailored specifically for any galaxy as would 
be the case during normal fitting. Vertical dashed lines are iso- 
metallicity lines with [Fe/H] = -1.3, -0.7, -0.4, 0.0 an 0.2 and 
horizontal solid lines are constant age 1.7, 2.8, 4.4, 7.0, 11.2, 14.1 
Gyr. 



Table 5. Derived stellar population parameters and stellar M/L 
ratios for the Coma compact galaxies with sufficient S/N. Upper 
error bounds set to max indicate the uncertainty estimate exceeds 
the models coverage in that direction. 



CcG 
(ID) 



Age (Gyr) 



[Fe/H] 
(dex) 



VI 

V9a 

V9b 

V19a 

V19b 



6.3±3.6 
12.0±2.1 

2.5±1.9 
11.9±2.4 
13.1±3.6 



p. -1 tr+macc 
^•-'^^-0.21 

-0.51±0.15 

n no+'max 
U.UC)_Q 25 

-0.12±0.1 
-0.56ib0.29 



[Mg/Fe 
(dex) 



(Mq/Lq) 



-0.04±0.16 
0.52±0.18 

0.3±0.3 
0.28±0.1 
0.48ib0.36 



5.6 
4.2 
1.95 
6.0 
4.2 



-\-max 

2.3 
+ 0.9 
0.8 
-\-max 
1.5 
1.2 
1.1 
+ 2.1 
1.2 



our sample. In all cases but CcG VI the default index weight- 
ing scheme of <Fe> and H/3 for [Fe/H] and age and Mgb 
for [Mg/Fe] was used. Other combinations were tested but 
results proved to be stable within errors. The fit of CcG VI 
failed with default index weighting because it marginally 
fell off the H/3-<Fe> model grid due to [Fe/H]> +0.2. As 
such we used only Fe5270 to constrain this galaxy's fidu- 
cial [Fe/H] and age. Uncertainties are estimated by repeated 
bootstrapping of the measured indices for each galaxy. We 
use the measured set of indices, perturbed them randomly 
by their respective error and ran them through the code. The 
la errors are then taken from the resulting parameter dis- 
tributions. This somewhat computer intensive approach was 
chosen, as opposed to using EZ-Ages' built in error scheme, 
because of situations where the task failed to compute uncer- 
tainties on the fiducial [Fe/H], age or both in either positive, 
negative or both directions. The issue apparently stems from 
the error bounds being off the model grid and can result in 
significantly under estimated uncertainties on for instance 
[Mg/Fe], whose derived error relies on those of the fiducial 
parameters and the Mgb index. Our method is still not ideal 
but provides a reasonable solution when trying to suitably 
assess error propagation. 
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Figure 11. Stellar population parameters for the Coma dwarf sample of S09 (blue points) and this work plotted against projected distance from 
the centre of the cluster (top row) and B band apparent magnitude (bottom row). The black dashed line represents the extrapolation of the 
relations between magnitude and stellar population parameters identified by S09. 



Also in Table [5] we present the inferred stellar M/L ra- 
ti os for our samp le. These are interpolated from Table 33 
of ISchiavonI (|2007l l and are based on scaled solar isochrones. 
The quoted errors are those introduced by the uncertainty 
on each galaxy's age and [Fe/H] and are also estimated by in- 
terpolation. The stellar M/L ratios do not take into account 
variable abundance patterns in the stellar atmospheres. 

Comparing the stellar and dynamical M/L for those 
galaxies in our sample with both, we see that the three cEs 
have dark matter fractio ns of ^ 40 to 70%. Th i s resu lt is 
similar to that found by IChilingarian MamonI (|2008l l for 
M59cO but, as one would expect, depends heavily on the 
choice of the virial coefficient [3. In addition, there is of 
course an aperture mismatch between that used to derive 
each of our M/L estimates in the sense that the Hectospec 
fibres are often smaller than 5Re. However, given that all the 
galaxies have colour profiles which either are flat or get bluer 
with radius, it is likely that the mean stellar M/L would 
at best remain constant and likely decrease if the aperture 
through which it was determined were enlarged. 

While it would now be desirable to compare our sam- 
ple's stellar parameters directly with data for other cEs from 
the literature, the inherent differences in the models used to 
derive such parameters makes this a difficult process. Ul- 
timately it is unclear whether the comparison would high- 
light differences in the models, stellar populations or both 
and disentangling them is be yond the scope of this work. 
However the recent work of ISmith et al. I (|2009l . hereafter 
S09) employed the same stellar population model as used 
here to analyse a sample of 79 faint red (g-r > 0.55) dwarf 
galaxies in the Coma cluster. Their results thus provide a 
unique comparative sample free of the systematic uncertain- 
ties mentioned above. 

In Fig. [11] we present the stellar population parameters 



of both data sets against projected radius from the centre of 
the cluster and B band magnitude. For this we take the cen- 
tre of the cluster to be RA=12:59:48, Dec=+27:58:50 and 
transform the SDSS magnitudes of S09 using the methods 
described previously. Considering the top row of panels first, 
in the [Fe/H] diagram the compacts are well mixed within 
the parameter coverage of the normal dwarf galaxies at their 
respective projected radius. The three bonafide cEs in our 
sample cover -0.6 ^[Fe/H]^ -0.1 with the two less extreme 
compacts, CcGVl and CcGV9b, having super solar metal- 
licity. In this plot, at least, the compact galaxies are largely 
indistinguishable from the sample of S09. The same cannot 
be said for the positions of the three clear cEs in the remain- 
ing panels which, despite the sizable errors, are all clustered 
towards older ages (> 12 Gyr) and higher [Mg/Fe] ratios (> 
0.25) than conventional dwarf galaxies of similar luminosity. 
This implies that the luminosity weighted stellar content of 
the cEs was formed rapidly at high redshift, a scenario which 
is generally associated with more massive galaxies. CcGVl 
and CcGV9b on the other hand have comparable ages to 
the mean of ^ 7=b3 Gyr seen for the S09 sample within 0.5 
Mpc of the cluster centre. 

Turning now to the lower row of panels of Fig. 1111 in 
the first diagram we see the relatively strong correlation, 
with slope —0.15 di 0.03, between magnitude and [Fe/H] for 
'normal' red Coma dwarf galaxies identified by S09. The 
following two panels present substantially weaker relations 
between magnitude and log(age) (slope 0.06 =b 0.03), and 
magnitude and [Mg/Fe] (slope —0.01 ±0.02). Unfortunately 
the magnitude overlap between this work's sample and that 
of S09 is rather limited making a robust comparison difficult. 
However, if we make the assumption that the slope of Fig. 
[2] is primar ily driven by metallic i ty, as has been concl uded 
elsewhere (jTerlevich et al.1 Il999l : IVazdekis et"all l200lh . we 
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can extrapolate the ms-[Fe/H] relation to fainter magni- 
tudes. It is then apparent that two of the Coma compacts, 
CcGVl and CcGV9b, have higher metallicities than would 
be expected even given the significant scatter in the relation. 
CcGV19a and CcGV19b are also located above the trend, 
but probably within the extension of the scatter, with only 
CcGVQa located on it. 



6 DISCUSSION 

As has been previously mentioned, the majority of theories 
surrounding cE formation rely on a deep potential tidally 
stripping some type of more massive progenitor galaxy and 
leaving only its more tightly bound, compact core intact. 
Logically the core should retain many of the properties of the 
central regions of the original galaxy, with some dependence 
perhaps on the severity of the tidal influence, and as such 
finding a compact object with these particular attributes 
gives us a direct clue as to its origin. With this in mind we 
review the findings of this work. 

(i) Six out of seven objects are found to have significantly 
redder B-Ic colours than the sequence defined by normal 
cluster members. Given how strong colour-magnitude corre- 
lations are known to be for cluster early- type galaxies, this 
clearly implies either a different formation mode to other 
cluster members or else an evolutionary difference. CcGV9a 
on the other hand is located within the scatter of the re- 
lationship. Crudely extrapolating back to the trend seen in 
Fig. [2]would indicate m^ ^ 18 mag or Mb ^ -17 mag, the 
top end of the dwarf regime, as the original magnitudes for 
the six outliers. Unfortunately given the fact that both E 
and early-type disk galaxies follow similar if not the same 
colour- magnitude relation in local galaxy clusters, it is im- 
possible to distinguish likely progenitors from this alone. 

(ii) Fitting PSF convolved point source + Sersic, single 
and double Sersic models we find that all of our sample 
are better represented by two components. Interestingly a 
similar result has also been found for a number of other 

20051 : 



cEs studied previously ( Graham 2002; Mieske et al, 

IChilingarian et al.1 lioOTl : Ichilingarian MamonI l2008l l and 



recently i n the best candidate f or a c E found in the Antlia 
cluster bv lSmith Castelh et al.] (|2008l l. 

Three of the sample, CcGVQa, CcGV19a and CcGV19b, 
are seen to occupy the same region on structural parameter 
diagrams as currently accepted cE galaxies (see Fig. [G] both 
panels) and so we conclude they should indeed be classified 
as such. Under the assumption cEs are of tidal origin, the 
expected direction of movement on the left panel in Fig.[6]of 
a Ms ~ -19 mag, </j.b >e'^ 20 magarcsec"^ and Re 1 kpc 
progenitor whose light profile is stripped from the outside in 
is likely close to the observed extension of the gE/bulge cor- 
relation. The right panel of Fig. [6] demonstrates that their 
inner components begin to fill in the parameter space be- 
tween Es/bulges and UCDs/nuclei. The other galaxy with a 
resolved inner component, CcGV9b, is less extreme in terms 
of surface brightness but is still substantially offset from the 
general trend seen in the right panel of Fig. [6] 

All of the remaining three galaxies are well fit by a point 
source + Sersic model with the absolute magnitudes of the 
point sources co mpa r able t o those found for galaxy nuclei in 
Virgo bv lCote et aD (|2006h . 



Fitting a gaussian to the nuclei 



luminosity function they obtain <Mg> = -10.74 mag (Mb ~ 
-10.4 mag) with a = 1.47 mag. Additionally they report a 
median half-light radius of 4.2 pc which would undoubtedly 
be unresolved at Coma, particularly given their magnitudes. 

(iii) Colour gradients are detected in four of the sample in 
the sense that their profiles become redder towards the cen- 
tre which, assuming there is little to no dust in these galax- 
ies, is evidence for radial variations in their stellar popula- 
tions parameters. CcGVl, CcGVlS and CcGV19b, are found 
to have little to no colour gradient but this does not rule out 
age or metallicity variations with radius. In general, signifi- 
cant changes in colour profile slope appear closely associated 
with the Re of the resolved cores of CcGV9a, CcGV19a and 
CcGV9b and an upturn is also noted for CcGV12. 

(iv) The velocity dispersions obtained for our sample are 
found to be greater than those of normal galaxies over the 
luminosity range in which the two populations overlap. 

Across —12 < Mb < —17 mag there is a notable corre- 
lation between magnitude and velocity dispersion for com- 
pact galaxies from both the literature and this work. Ex- 
cluding CcGV9b due to its significant deviation from the 
trend and considering ground based observations only, we 
obtain x = 3.7 when fitting L oc across this range. It 
would be difficult to assess the uncertainty on x due to the 
inhomogeneity of the data but the similarity to that found 
previously for Es and bulges, x ^ 4 (Faber Jackson 19 76), 
as opposed to dEs, x ^ 2 (|Matkovic GuzmanI 120051 1 , is 
noteworthy. Additionally, the offset nature of this relation 
is perhaps one of the strongest pieces of evidence for the 
progenitors of cEs being more massive galaxies which have 
undergone some form of stellar mass removal. It is worth 
pointing out that this would imply taking an early-type 
galaxy, using some method to dim it by 1-3 magnitudes (or 
equivalently, assuming a constant stellar M/L, removing be- 
tween 60% and 95% of its stellar mass) and it s till re main 
intact. This issue was highlighted bv lBekki et al.1 (|200lh who 
performed hydrodynamical simulations of a satellite galaxy 
interacting with M31. They concluded that the scarcity of 
cEs can be put down to the small range of initial orbital 
parameters that allow, in their study, a spiral galaxy to be 
transformed via galaxy "threshing" into a cE and that once 
formed the lifetimes of cEs would be short (less than a few 
Gyr) before being full accreted. Their results also estimate a 
^ 60% mass loss for the disk component of the satellite but 
only a ~ 20% mass loss for the bulge. Such a reduction in 
the bulge stellar content would be insufficient to produce the 
observed position of many cEs seen in Fig. [9l Another pre- 
diction of their model is the triggering of a nuclear starburst 
in the cE remenant of the satellite galaxy. The three cEs in 
Coma show no evidence for this based on their colour pro- 
file slopes and particularly mean luminosity weighted ages 
within the Hectospec fibres. 

Alternatively, Moore et all (jl996| ) report that cEs are a 
rare class of elliptical galaxy that can form alongside normal 
dEs via the galaxy-galaxy harassment of bulge-less spiral 
galaxies. However, explaining the two component nature of 
many cEs would seem difficult in this scenario. 

As a third possibility, some cE galaxies could be extreme 
examples of "hyper-compact stellar systems" (HCSSs), de- 
fined as clouds of stars bound to supermassive black holes 
that have been ejected from the centres of their host galax- 
ies (Merritt, Schnittman and Komossa 2008). Effective radii 
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and measured velocity dispersions of HCSSs are predicted 
to be 

^ 50pcf-^V ^^V', (5a) 

^ \10^Mq J V500 km s-V 

cr ^ 100 km s f ^ T ) (5b) 

V500 kms-V 

where Vk is the kick velocity and M, is the black hole mass. 
Explaining the Coma cE galaxies in terms of HCSSs would 
require large (> 10^ M©) black holes and kicks below galactic 
escape velocities (Vfc < 500 km s~^). 

(v) Employing up to date SSP models and fitting code 
we have constrained the luminosity weighted age, metallic- 
ity and [Mg/Fe] ratio of our sample within the ^ 700 pc 
diameter of the Hectospec fibres at Coma. The three cEs all 
have old ages (> 12 Gyr), intermediate metallicities (-0.6 < 
[Fe/H] < -0.1) and high [Mg/Fe] (> 0.25). Such parameters 
are often associated with more massive galaxies and particu- 
larly places any in-situ formation scenario, that is to say not 
the transformation of some progenitor, in doubt given the 
expected lifetimes of cEs. Of the two other Coma compacts 
with sufficient S/N spectra, both have young to intermediate 
ages and solar or super-solar metallicities albeit with large 
error bounds. CcGVl is found to have solar [Mg/Fe] while 
CcGV9b has [Mg/Fe] = 0.3 which is unusually high for its 
age, although again there is the caveat of a large uncertainty. 

We find that four out of five galaxies have [Fe/H] either 
higher (CcGVl and 9b) or at the upper limit (CcGVlQa and 
19b) of that expected for their luminosity. As such we can 
again extrapolate back and note that a 1-3 magnitude shift 
would be required to bring them back in line with the bulk 
of normal galaxies. CcGV9a's position on the extreme of the 
trend in Fig. [2] is likely generated by a combination of its 
old age and metallicity characterised by low [Fe/H] and high 
[Mg/Fe]. 

While somewhat unwise to directly compare stellar pop- 
ulation parameters derived from different models, it is at 
least interesting to note those associated with the previ- 
ously known cE population. A recent study by iRose et al] 
(|2005l ) found M32 to have an age and metallicity variation 
from 4 Gyr and [Fe/H ] — 0.0 in the nucleus to 7 Gyr and 
[Fe/H] = -0.25 at li^.. I Sanchez- Blazauez et al.1 (l2006bl ) find 
NGC 4486B to have age 9.5 Gyr and [Fe/H] = 0.4 and 
NGC 5846A to have age 10.7 Gyr and [Fe/H] = 0.08. Fi- 
nally, [Ctiilingarian et al.l (|2007f ) and IChilingarian MamonI 
(|2008l ) obtained ages and metallicities for A496cE and 
M59cO of 16.4 Gyr and [Fe/H] = -0.04 and 9.3 Gyr and 
—0.03 respectively. Essentially, the trend appears to be to- 
ward intermediate to high ages and metallicities which is in 
agreement to the parameters found for our cEs. In terms of 
[Mg/Fe], A496cE and M59cO have super solar values (-- 0.2) 
while M32 has sub solar, likely a result of its relatively young 
luminosity weighted age. 

Of course one major issue here, aside from model depen- 
dence, is differences in instrumental setup and aperture size 
given the strong gradients found in M32, at least. Also M32 
has significantly better spatial resolution than the other 
galaxies which m akes suitable compar isons difficult. How- 
ever, we note that lPavidge et all (|2008l ) conclude from their 
high resolution K band spectroscopic observations of M32, 
NGC 4486B and NGC 5846A that the former has the chem- 



ical maturity of a spiral galaxy and is less Fe deficient than 
the latter pair. 

(vi) For those galaxies in our sample which have both a 
measured dynamical and stellar M/L we estimate that they 
have dark matter mass fractions ~ 40% to 90%. Our colour 
profiles indicate this result should be robust to the aperture 
mismatch between the Hectospec fibre (radius ^ 1 to 5i?e) 
and the radius used to derive our dynamical mass estimates 
(5i?e)- However, particularly in the case of CcGVl, the de- 
gree to which this figure may be effected by the inherent 
uncertainties in deriving dynamical M/L for stellar systems 
which are probably not virialised is unclear. 



6.1 Progenitors 

Combining the above conclusions, the present work sup- 
ports, to varying degrees, the tidal disruption and strip- 
ping scenario as the formation mode of compact galaxies, 
CcGV12 and CcGV9b being the weakest contributors owing 
to their poorly constrained velocity dispersions. Assuming 
this is the correct model, progenitor morphology is never- 
theless difficult to assess. 



6.1.1 CcGVOa 

Of the clear cEs, CcGV9a seems to have at least two defining 
similarities to those of typical disk galaxies in the presence of 
a weak bar and possible outer ring. The outer component's 
light profile, with n — 1.64, might perhaps then represent 
the threshed remains of its disk. As such this galaxy is our 
best candida te to fit a similar f ormation mechanism to that 
suggested bv lBekki et al.l (|200l[ ) and have an early- type disk 
galaxy as its progenitor. Furthermore, based on the signifi- 
cant fraction of light contained within its outer component 
and relatively strong signal of its original morphology, we 
speculate that this galaxy may still be in the early stages of 
being transformed into a cE. On the other hand, how such a 
strong colour gradient is generated over such a small spatial 
scale in a class of galaxy which conventionally has rather 
shallow gradients is unclear. In fact the observed colour gra- 
dient is in the opposite sense to that predicted by the model, 
and a luminosity weighted age of ^ 12 Gyr and [Fe/H] = 
—0.5 within the Hectospec fibre largely rule out a stripping 
induced starburst in the centre. 



6.1.2 CcGVlQa 

The likely progenitor of CcGVl 9a is even less well defined. It 
has a number of properties in common with CcGV9a yet an 
inner Re — 76 pc would seem rather small for a normal bulge 
and would have to be the result of it shrinking during the 
tidal stripping phase. Alternatively, the inner component of 
this galaxy is only just larger than the upper limit found for 
elliptical galaxy n uclei [Re — 62 pc) in the Virgo cluster by 
I Cote et al.l (|2006l V Ultimately there is no strong pointer ei- 
ther way except perhaps the assumption that SO-Sa galaxies 
may possess stronger colour gradients than elliptical galax- 
ies. 
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6.L3 CcGVlOb 

Our best assessment for the final definite cE, CcGV19b, 
would be that it is the stripped remains of a -17 mag 

dE,N given its compact inner component, flat colour pro- 
file, velocity dispersion and stellar population parameters. 
Furthermore, this galaxy has properties which are close, al- 
beit more luminous and more massive, to those of UCDs 
in For nax and Virgo (iMieske et al.ll2006l:[ Evstigneeva et al.l 



l2007bl l. e.g. UCD 3 in Fornax ('Phillip DS et al.1 [20011 ) wh'ich 

has a compact core and relatively large outer halo, and so 
we speculate this galaxy may be an example of a UCD pro- 
genitor. 

6.14 CcGVl, CcGV12, and CcGVlS 

All three galaxies from our sample which have unresolved 
nuclei and outer components with n — 1-1.9, more than 
likely originated from nucleated elliptical galaxies with ini- 
tial yis'^ -17 mag. In addition, CcGVl certainly has stellar 
population parameters which indicate substantial gas pro- 
cessing and only relatively recent star formation truncation 
which in turn may imply recent entry to the core of the 
cluster (see S09). 

While strictly speaking not sufficiently extreme in terms 
of surface brightness to be classified as cEs or rather M32- 
like, these three galaxies are all more spatially compact than 
conventional dwarf galaxies at their luminosity and have 
higher velocity dispersions. It would appear likely they orig- 
inated as the result of either less severe tidal disruption, 
different progenitor properties, or both, relative to cEs. 

6.1.5 GcGVOb 

CcGV9b is an odd ball even within this sample. Extrapolat- 
ing back the position of its inner component on Fig.[6]points 
to the dE region of the plot, assuming it experiences similar 
scenarios to the rest of the sample, and it has a velocity dis- 
persion, albeit poorly constrained, which compliments this 
conclusion. Yet its inner component is also seemingly too 
large. Re — 145 pc, to be a classic nucleus and the galaxy's 
colour profile is strongly sloped. Additionally, this galaxy 
has the youngest luminosity weighted age of the sample. At 
best we can comment that this galaxy must have only re- 
cently been introduced to the inner regions of the cluster 
but cannot robustly draw any further conclusions from our 
data alone. 



7 SUMMARY 

In this paper we have reported the discovery of seven com- 
pact galaxies in the Coma cluster that begin to fill in the pa- 
rameter space between UCDs and cEs. Three were initially 
identified from HST/ACS imaging and MMT/Hectospec 
spectroscopy with a further four candidates being eyeball 
classified as high probability cluster members with support 
from their colours, size and resolved appearance in the ACS 
imaging. Subsequent spectroscopy confirmed all four candi- 
dates to be cluster members. 

All seven compact galaxies are found to be well fit by a 
two component light profile and have structural parameters 



that significantly deviate from those expected for galaxies at 
their luminosity or size. The measured velocity dispersions 
for our sample exceed those of conventional dEs galaxies and 
the three extreme enough to be classified as cEs are found to 
have old, intermediate metallicity stellar populations with a 
substantial Fe deficiency relative to solar abundance. 

Our best explanation is that the compact galaxies in 
this work originated as more massive systems that were sub- 
ject to some form of tidal mass removal and that they, partic- 
ularly the cEs, represent the rare, extreme outcome of such 
galaxy-galaxy processes. We also speculate that the progeni- 
tors of our sample come from at least two distinct early-type 
morphologies, both of which have suffered similar nurturing 
by the dense cluster environment. 

In Ught of the success reported here which employed 
relatively simple criteria, we propose that further searches 
should be mounted in other massive clusters to continue to 
study this extreme end of galaxy formation and evolution. 
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